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The origin of cosmic ray events with E > 10 11 GeV remains mysterious. In this 
talk I briefly summarize several proposed particle physics explanations: a breakdown of 
J> \ Lorentz invariance, the "Z— burst" scenario, new hadrons with masses of several GeV 

as primaries, and magnetic monopoles with mass below 10 10 GeV as primaries. I then 
q . describe in a little more detail the idea that these events are due to the decays of very 

massive, long-lived exotic particles. 
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1 Introduction 



The observation [1] of cosmic ray (CR) events with energy E > 10 11 GeV (ultra-high energy, 
UHE) poses at least two distinct problems: 

• The energy problem: It is very difficult, if not impossible, to accelerate protons to 
such energies in any known astrophysical source [2]. 

• The propagation problem: Protons with E > 5- 10 10 GeV can photoproduce pions on 
the cosmic microwave background (CMB). The GZK effect [3] implies that protons with 
E > 10 11 GeV should have been produced within a few dozen Mpc of Earth. The same 
conclusion holds for photons [which get absorbed via 7uhe + 7radio — > e + e~, where 7 ra dio 
belongs to the (extra)galactic radio background] and heavier nuclei (which are broken 
up by collisions with CMB photons). Current estimates of (extra)galactic magnetic field 
strengths imply that protons with E > 10 10 GeV should point back to their source if 
they are produced within a few dozen Mpc. There are, however, no known near sources 
of high-energy particles in the direction of the most energetic events.^ 

Many suggested solutions of "the UHECR puzzle" actually only address the second problem. 

1 will briefly review these ideas in the following section. In sec. 3 I will discuss two approaches 
that can solve both problems, the main emphasis being on decays of super-massive particles, 
before concluding in sec. 4. 

2 Partial solutions 

2.1 Breaking Lorentz invar iance 

Calculations of the GZK effect assume that u^p — > Nir) (N = n,p) is the same for E 1 ~ 
Equb ~ 10" 3 eV, E p ~ 10 11 GeV and E 1 ~ 200 MeV, E p = m p . This assumption may 
be wrong if Lorentz symmetry is violated. This is actually quite an old idea [6], which has 
been rediscovered recently [7] . By assuming different limiting velocities for different species of 
particles one can arrange for a suppression of the GZK effect; one can even avoid it altogether. 
In this case one would expect a completely smooth spectrum around the GZK cut-off. In 
particular, one would not expect a bump just below the GZK energy, which seems to be 
present in both the AGASA and HiRes results [1]. Besides, it seems to me a step backwards 
to give up one of the basic symmetries underlying our understanding of Nature; it rather 
reminds me of Bohr's willingness to give up energy conservation to explain nuclear f3 decay 
spectra. Back then, Pauli's neutrinos saved the day; I'm confident that the puzzle of the 
UHECR events can also be solved without giving up well-established principles. 

2.2 Exotic hadrons as UHECR primaries 

If an exotic hadron He with mass rrih E = rm p is [8] the primary of the UHECR events, the 
GZK "cut-off" (better: spectral break) will be pushed upwards by a factor r. Even the 
AGASA data can be explained if r > 4. On the other hand, if Yie is too heavy, it looses 

'''There have been claims [4] of statistically significant correlations between the arrival directions of UHECR 
events and certain kinds of active galactic nuclei (AGN), thought to be possible sources of UHE particles. 
However, an independent analysis did not confirm this [5]. 
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energy too slowly in the air shower it initiates when hitting the Earth's atmosphere. It has 
been estimated [9] that the produced shower will only look sufficiently proton-like if r < 50. 

There are at least two problems with this scenario. First, acceleration becomes much 
more difficult to explain. Producing UHE He particles by first accelerating protons and then 
colliding them with ambient matter (or photons) will likely over-produce neutrinos and UHE 
photons, since the (inclusive) cross section for pion production is several orders of magnitude 
larger than the cross section for He production; of course, this would also assume that protons 
can be accelerated to energies well beyond the most energetic observed CR event. Directly 
accelerating stable Ke particles is also problematic, since at least on Earth they only constitute 
at most a tiny fraction of all matter, so the heavenly accelerator would presumably also "waste" 
most of its energy on accelerating ordinary protons and electrons. The second problem is 
that He particles, being strongly interacting and rather light, should have been produced 
abundantly in collider experiments; it is difficult to believe that they could have escaped 
detection. 

2.3 Z bursts 

The idea [10] is that UHE neutrinos are produced somewhere; since they can propagate almost 
freely through the Universe, this neutrino source may be anywhere within one Hubble radius, 
3 Gpc, or so. Within the GZK radius around the Earth, a small fraction of these UHE 
neutrinos annihilate on relic background antineutrinos to produce Z bosons, most of which 
decay hadronically, giving rise to the observed UHECR events; of course, UHE v annihilating 
on relic v also contribute. 

One problem with this explanation is that it aggravates the energy problem by several 
orders of magnitude. One will need UHE neutrinos of at least five times the energy of the 
most energetic UHECR event, since Z bosons decay into typically > 20 hadrons. In order to 
produce such neutrinos, one needs protons whose energy is higher by at least another factor 
of five or so, i.e. one would need a source of protons extending to E > 10 13 GeV. Moreover, 
the flux ^uhe 3> ^exp. uhecr; i- e - the intensity of this source must be much higher than that 
needed to directly accelerate the observed UHECR flux.* Indeed, the required UHE neutrino 
flux is at best marginally compatible with existing limits [11]. 

3 Complete solutions 

3.1 Magnetic monopoles as UHECR primaries 

This is another old idea [12] that has been rediscovered recently [13]. The main observation 
is that magnetic monopoles can actually be accelerated (rather than only deflected) by (ex- 
tra)galactic magnetic fields. Indeed, this acceleration is quite efficient, partly due to the large 
effective charge (> l/a cm ) of the monopoles. According to current estimates energies beyond 
10 11 GeV are easily accommodated in this way. Of course, monopoles are stable, so they can 
have been produced in the very early Universe, most plausibly during a phase transition. The 
density, or flux, of monopoles required to explain the observed UHECR events is safely below 
all known bounds [14], as long as these monopoles do not catalyze nucleon decay. 

*Such a source could easily over-produce protons just below the GZK cut-off. 
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The biggest difficulty of this model is to explain why the observed events look so much 
like proton-induced air showers. In order to produce a shower at all, the monopole must be 
ultra-relativistic, i.e. m M < 10 10 GeV. These monopoles can therefore not be associated with a 
Grand Unified symmetry; instead one has to postulate the existence of an "intermediate" scale 
Oirau)- Of course, direct monopole searches at colliders imply tjim 3> m p . Simple kinematics 
then implies that a monopole-induced air shower will penetrate much more deeply into the 
atmosphere than a p— initiated shower does, unless the cross section for inelastic scattering of 
a monopole on an air nucleus is much larger than the corresponding cross section for protons. 
The authors of ref.[14] argue that bound states of monopoles carrying nontrivial SU(3) C charge 
might conceivably have this property. In this view, successive collisions could increase the 
length of the color-magnetic QCD string (or flux tube) between the constituents, and hence 
the area of the bound state. Note that this string can only break into monopole-antimonopole 
pairs. Since monopoles are heavy, a purely color-magnetic string might therefore store many 
orders of magnitude more energy, and could correspondingly be many orders of magnitude 
longer, than the more familiar color-electric string does, which can break up into light qq 
pairs. The effective monopole-air cross section can grow quickly through this mechanism, 
leading to the desired short air shower. Note, however, that the monopoles, being stable, 
should still reach the Earth (with non-relativistic velocity); detecting these monopoles would 
be the ultimate test of this explanation. 

3.2 Decaying superheavy particles 

Another exotic possibility is that the UHECR events are due to the decay (or annihilation) 
of some super-heavy long-lived X particle [15]. We obviously need m x > 10 12 GeV in order 
to generate a significant flux of UHECR primaries at E > 3 • 10 11 GeV. Just as obviously, 
the lifetime of X must be at least comparable to the age of the Universe, Tx ^ 10 10 yrs. 
This immediately raises two particle physics questions: how can such massive particles be 
so long-lived, with T x /m x < 1CT 54 ? And how could such massive particles ever have been 
created? 

One way to ensure the longevity of the X particles is to embed them into topological defects. 
The original proposal [15] envisioned annihilating monopole-antimonopole pairs as source of 
UHECR events. This idea has later been refined to the concept of a "cosmic necklace", 
( ant i) monopoles strung along a cosmic string [16]. Cosmic strings themselves become unstable 
when they intersect. Such kinds of defects could be formed during a phase transition via the 
Kibble mechanism [17]. 

Note that we now need the scale of symmetry breaking to be > 10 12 GeV (as opposed 
to < 10 10 GeV in the previous subsection). This may be problematic, since with the scale 
of (four-dimensional) inflation determined to be ~ 10 13 GeV, it is quite possible that the 
post-inflationary Universe never was hot enough to have been in the "unbroken" phase, which 
would require T > 10 12 GeV. Moreover, the increasingly well measured CMB anisotropics do 
not seem to show any evidence that topological defects (like cosmic strings) played a role in 
structure formation. In these topological defect models the X particles could be superheavy 
gauge or Higgs bosons, and/or their superpartners. 

Alternatively the X particles could exist as (more or less) free particles in today's Universe. 
Such particles might have been created at the end of inflation, when the energy density of 
the Universe was more than one hundred orders of magnitude higher than it is today. Several 
production mechanisms have been suggested, e.g. gravitational production due to the rapidly 
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varying metric (which could work for X particles as heavy as 10 16 GeV) [18], or in (inclusive) 
inflaton decays [19]. The X particles would then serve as "batteries", storing energy from this 
extremely violent early epoch of the Universe and releasing it in our much balmier times. 

In order to explain the longevity of such free X particles one has to postulate that their 
couplings to ordinary matter are very strongly suppressed. Some such suppression is actually 
expected if X resides in the "hidden sector" thought to be responsible for the spontaneous 
breaking of supersymmetry [20], which by definition only has M P i anck suppressed couplings to 
the visible sector, although it is by no means guaranteed that this suppression is sufficient. The 
relevant couplings could also be suppressed by (approximate) symmetries [21], or geometrically 
in brane world scenarios [22].* 




1 TeV 1 GeV 

(SUSY (hadronization) 
+ SU(2) ® U(l) 
breaking) 



Figure 1: Schematic MSSM cascade for an initial squark with a virtuality Q ~ Mx- The full circles 
indicate decays of massive particles, in distinction to fragmentation vertices. The two vertical dashed 
lines separate different epochs of the evolution of the cascade: at virtuality Q > Msusy, all MSSM 
particles can be produced in fragmentation processes. Particles with mass of order Mgusy decay at 
the first vertical line. For Msusy > Q > Qhad light QCD degrees of freedom still contribute to the 
perturbative evolution of the cascade. At the second vertical line, all partons hadronize, and unstable 
hadrons and leptons decay. See the text for further details. 

The spectrum of UHECR primaries (and other stable particles) at source is determined by 

*It is sometimes claimed [23] that the required long lifetime of X needs severe finetuning. This is not really 
true, since tx could be as large as 10 20 yrs [24] if X particles constitute the bulk of Dark Matter. Model 
builders therefore have some ten orders of magnitude in tx to aim for. Of course, the observed UHECR flux 
determines the ratio fix / tx to within a factor of a few, but this does not require any more finetuning of model 
parameters than in any other explanation of the UHECR events. 
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the physics of X decays [25]. This is quite nontrivial, as indicated in Fig. 1, which is taken 
from ref.[26]. The primary few-body X decay initiates a generalized parton shower. This 
is similar to, e.g., the hadronic decay of Z bosons, which have been studied in great detail 
at LEP and SLC. From a theorist's perspective, a Z boson decays into a qq pair, whereas 
experimentally one observes twenty or so hadrons in the final state. The transition from the 
two (anti) quarks to the O(20) hadrons entails both perturbative and non-perturbative physics. 
The perturbative part is called a parton shower. It can be envisioned by assuming that the 
original qq pair is not produced on-shell, but with initial time-like virtualities of order m Z - 
The showering occurs when these partons move closer to the mass shell by emitting additional 
partons, mostly gluons, with smaller time-like virtualities. 

Similarly, the particles produced in primary X decay should be considered to have initial 
time-like virtualities O(mj). At such high energy scales parton showers are expected to de- 
velop quite differently from parton showers at scale m Z - To begin with, the existence of a (real, 
four-dimensional) energy scale mx rnz strongly indicates the existence of superparticles, 
which can stabilize this very large ratio of scales against radiative corrections [27]. Since the 
superparticle mass scale msusY i$ 1 TeV <C mx, superparticles can be produced (i.e., will be 
"active") in the initial part of the shower development, even if X decays primarily into SM 
particles. Moreover, we know that at energy scales near m x all three gauge interactions are 
of approximately equal strength. In refs.[26, 28] we therefore included all gauge interactions 
as well as third generation Yukawa interactions in the description of the parton shower. 

When the shower virtuality scale reaches msusY, superparticles decouple from the shower 
and decay; the electroweak gauge and Higgs bosons, as well as top quarks, decouple at es- 
sentially the same scale. A careful treatment of these decays is mandatory if one wants to 
account for the entire energy released in X decays. As well known, sparticle decays can be 
somewhat complicated even at the parton level, involving lengthy decay cascades. We treated 
these with the help of ISASUSY [29]. 

At shower scales between msusY and Qhad ~ 1 GeV only standard QCD is important for 
the evolution of the parton shower, as at LEP. At scale Qhad the non-perturbative parton — > 
hadron transition occurs. Most hadrons, as well as heavy \x and r leptons, eventually decay 
as well, as indicated in the figure. At the end one is thus left with seven species of particles 
(plus their antiparticles): protons, electrons, photons, three kinds of neutrinos, and LSPs. 

Out of the multitude of particles produced in a far-away X decay at most one will be 
observed on Earth. We therefore "only" have to compute the one-particle inclusive X decay 
spectra into the seven stable decay products. In QCD the notion of fragmentation functions 
(FFs) has been developed precisely in order to describe single-particle inclusive spectra [30]. 
These FFs depend on two variables: the energy scale of the shower, here m^, and the scaled 
energy xp = lEpjmx of the stable decay product P in the X rest frame. The x dependence of 
hadronic FFs at some reference scale has to be taken from data, since non-perturbative effects 
are important here; we used the fits of ref.[31].* The scale dependence of the FFs is described 
by the appropriate generalization of the well-known DGLAP evolution equations [32]. The 
decays of superparticles and other massive partons can also be treated in this framework, as 
long as their energy is much larger than their mass, in which case a collinear treatment is 
adequate. Finally, at very small x color coherence effects should be included [26]. 

The main results of this analysis can be summarized as follows [28, 26]: 

*Some of these "input" FFs are still not well determined, e.g. those starting from a b or c quark. 



6 



• For small x < 0.01, the z/ e , z/ M , e, 7 and p fluxes^" have essentially fixed ratios, independent 
of the primary X decay mode(s) and of the details of the SUSY spectrum. The flux 
is largest, exceeding the p flux by a factor ~ 3, while the photon and e ~ u e fluxes 
exceed the proton flux by a factor ~ 2.5 and ~ 2, respectively. All these fluxes can 
be described by a simple power law in this region, d^/dE oc E~ lA . These features are 
due to QCD evolution effects. Note that most z/ ejM , e and 7 at x < 0.01 originate from 
pion decays, and are thus subject to the same QCD effects. * On the other hand, the v T 
and lightest superparticle (LSP) fluxes, which decouple much earlier in the shower, are 
model-dependent even at small x. Generally they increase more slowly with decreasing 
x, and thus become subdominant for x < 0.01. 

• For x > 0.1 the fluxes at source depend very strongly on the primary X decay mode(s), 
and less strongly (typically at the factor of two level) on details of the SUSY spectrum, 
the relative ordering of and mass splitting between states being more important than 
the overall sparticle mass scale. For all X decays into MSSM particles the photon flux 
at source is larger than the proton flux. In case of hadronic X decays the ratio <& 1 /<& p 
reaches a minimum (slightly above 1) at x ~ 0.3. At larger x it increases due to the 
emission of hard photons early in the shower, illustrating the importance of including 
electroweak interactions in the description of the shower. At smaller x this ratio increases 
because the ratio of FFs into pions and into protons increases with decreasing x [31]. 
For primary X decays into only weakly interacting particles the photon flux at x > 0.1 
exceeds the proton flux by at least one order of magnitude, since in this case hadrons 
can only be produced after two electroweak branching processes, while photons can be 
emitted already in the first branching. The and v e fluxes also exceed the proton flux 
at large x; these two neutrino fluxes approach each other at x>0.5, where most neutrinos 
come directly from electroweak branching and decay processes, rather than from pion 
decays. 

• The LSP flux is important even if the primary X decay only involves SM particles. For 
example, Bino-like LSPs carry 2-3% of the total energy if X decays into leptons, 5-6% 
if X decays into quarks, 20-30% if X decays into squarks, and 30-50% if X decays into 
sleptons. The percentage is higher for X decays into SU(2) singlet sfermions, since they 
have shorter supersymmetric decay chains, and also shower less (due to the absence of 
SU(2) couplings). 

These fluxes are modified by propagation effects. Even if most relevant X decays occur 
in the halo of our own galaxy, which is expected [33] if X particles can move freely under 
the influence of gravity, electrons will loose most of their energy in synchrotron radiation on 
galactic B— fields. Moreover, the propagation distance is so large that oscillations essentially 
wash out all information about the original neutrino flavor dependence, i.e. approximately 
equal fluxes of v e , and v T will reach Earth independent of the ratios of these fluxes at source. 

How can one test this class of models? Given that m x and the source density oc Qx/ r x 
are considered to be free parameters, it is not surprising that one can fit [25, 34] the observed 
UHECR flux. These models are even compatible with a GZK spectral break, as indicated by 

t All fluxes are summed over particle and antiparticle, if they are distinct. 

*These effects also imply that the results are largely independent of the necessary extrapolation of the 
"input" FFs towards small x, as long as this extrapolation conserves energy and has a milder small— x behavior 
than x~ 1A . 
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the current HiRes (but not AGASA) data, if X particles are distributed more or less uniformly 
throughout the Universe, as would e.g. be expected if they are confined to cosmic strings. 
As mentioned earlier, freely moving X particles are expected [33] to have a large overdensity 
(~ 10 4 to 10 5 ) in our galaxy, in which case most relevant X decays would occur at distance 
well below one GZK interaction length. Even in this case the HiRes spectrum can be described 
by X decays, simply by choosing mx such that the spectrum cuts off just above 10 11 GeV. 

The composition of UHECR primaries is more problematic, however, Observations indicate 
that most primaries are protons or heavier nuclei, not photons. This follows from an analysis 
[35] of the longitudinal shower profile of the most energetic event ever observed; from the 
large number of muons in the Haverah Park [36] and AGASA [37] data; and from the absence 
of a South-North asymmetry in the AGASA data [37], which would be expected for photon 
primaries, since photons with E 1 > 10 10 GeV split into e + e~ pairs already high in the Earth's 
magnetosphere. In contrast, as mentioned earlier, X decay models predict the photon flux 
at source to be higher than the proton flux. If most sources are at cosmological distances, 
propagation effects may well suppress $ 7 at post-GZK energies even more than $ p . However, 
if most sources reside in the halo of our galaxy, propagation effects on $ 7 are expected to be 
significant only if the galactic radio background has been underestimated by at least an order 
of magnitude. I do not know how (im)plausible this might be. 

Fitting the observed UHECR flux to $ p alone increases the predicted flux of neutrinos 
and LSPs on Earth by about a factor of three if m x 3> 10 12 GeV; if m x ~ 10 12 GeV, the 
enhancement amounts to a factor ~ 2 for hadronic primary X decays, and to roughly one order 
of magnitude for leptonic decays. § Even without this enhancement the neutrino flux might 
be detectable by future experiments [38] . One can look for long muon tracks as well as high- 
energy particle showers due to NC and u eiT CC events in neutrino telescopes like IceCube [39] . 
Another technique, used by the RICE experiment at the South Pole [40], is to look for coherent 
Cherenkov emission of radio waves, which is mostly sensitive to electromagnetic showers, i.e. 
to u e CC events. Finally, future air shower detectors like the Pierre Auger array [41] can look 
for "horizontal" air showers, which originate too deep in the atmosphere to be due to photons 
or protons, and for r— leptons coming out of the Earth at shallow angles and decaying in the 
atmosphere. In all cases one has to require E > 100 TeV in order to suppress the atmospheric 
neutrino background. 

Some estimated [34] event rates per year are collected in Table 1, for two extreme values 
of mx and two extreme X distributions. "Galactic" here means that essentially all relevant 
sources are closer than one GZK interaction length. In contrast, if X particles are distributed 
homogeneously throughout the Universe, X decays occuring at distances well beyond one GZK 
interaction length contribute to $„, but not to $ p (at post-GZK energies). A homogeneous 
distribution therefore leads to ~ 10 times higher neutrino event rates than a galactic distri- 
bution does. Similarly, mx ~ 10 12 GeV leads to an about ten times higher neutrino event 
rate than mx ~ 10 16 GeV does. The reason is that, in particular for a galactic distribution 
of X particles, models with mx ^> 10 12 GeV only allow to fit the very highest energy end 
of the UHECR spectrum, leading to significantly smaller fluxes at source. Notice also that 
the cross sections for up scattering grow significantly less rapidly than linearly with energy 
once E v 3> 1 TeV [42]. Finally, the ranges in the table indicate the spread between different 
X decay models. For m x ~ 10 16 GeV all relevant particles are at a; < 1, in which case this 

§Leptonically decaying X particles with m x < 10 13 GeV are thus particularly difficult to reconcile with the 
experimental evidence against photons as main UHECR primaries. 
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model dependence disappears, as mentioned above. 



Table 1: Number of neutrino events with E v > 100 TeV per year expected in different detectors. 
See the text for further details. 



m x [GeV] 


X distribution 


IceCube 


RICE 


Auger 


2 • 10 12 
2 • 10 12 


Galactic 
Homogeneous 


10-30 
80-300 


1-4 
10-35 


1-3 
10-25 


2 • 10 lb 
2 • 10 16 


Galactic 
Homogeneous 


1 

10-15 


0.4 
6 


0.3 
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Unfortunately the prediction of a detectable neutrino flux at energies beyond that reached 
by atmospheric neutrinos is hardly unique for this explanation of the UHECR events. The 
very GZK effect itself gives rise to extremely energetic neutrinos through the decay of charged 
pions. Note also that the numbers in table 1 differ by more than two orders of magnitude, i.e. 
a large range of neutrino fluxes is broadly compatible with X decay models. These models 
generically predict that the neutrino flux extends to somewhat higher energies than the proton 
flux does [26]. However, testing this prediction would require a fairly accurate measurement 
of the neutrino flux at energies well beyond 10 11 GeV, where the expected event rate in the 
experiments covered in table 1 is very low. 

An unambiguous test of this kind of model could be achieved by detecting the flux of very 
energetic LSPs. In "bottom-up" scenarios, where the observed UHECR events are explained 
by the acceleration of protons (or heavier nuclei), very energetic superparticles could only be 
produced when the accelerated protons (or nuclei) collide with ambient matter. However, the 
inclusive cross section for pion production is at least seven orders of magnitude larger than 
that for the production of superparticles. In this kind of model a detectable LSP flux would 
therefore be accompanied by a huge flux of extremely energetic neutrinos, which would be 
(comparatively) easily detectable. In contrast, we saw earlier that LSPs are expected to carry 
a significant fraction of the energy released in X decays. 

Detecting this UHE LSP flux is difficult, but may not be entirely hopeless [43]. The 
cross section for LSP-nucleon scattering is estimated [44] to be 10 to 100 times smaller than 
the neutrino-nucleus cross section, if the LSP is Bino-like, which is the case in most SUSY 
models. These LSPs will therefore be absorbed less strongly in the Earth than neutrinos are. 
For E v = 10 9 GeV only one in a thousand neutrinos impinging on Earth at an angle at least 
5° below the horizon can transverse it without interacting. u T s can be regenerated via r — > v T 
decays even if they undergo CC scattering, but the emerging v T will only carry 20% or so of the 
energy of the original one. If the UHE neutrino "background" comes from X decays, one can 
therefore show that very few events with E > 10 9 GeV can emerge at an angle 9 > 5° below 
the horizon, and essentially none at 9 > 10°. In contrast, if a^ p = a up /10, at E% = 10 9 GeV 
the LSP flux will be depleted significantly only at 9 > 60°; the depletion becomes altogether 
negligible for o^ p = a up /100, up to energies well beyond 10 10 GeV. 

Due to the LSPs' small flux and cross section, one will need huge targets. These might be 
provided by space-borne fluorescence detectors like EUSO [45] and OWL [46]. Table 2 shows 
estimated [43] rates of LSP events with 9 > 5° and E% > 10 9 GeV, assuming a^ p = a up /10 
and a target area of 150,000 km 2 . Since LSPs can interact either in the atmosphere or in ~ 10 
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m water-equivalent (w.e.) below it, this corresponds to an approximate target size of 2,000 
km 3 w.e. . Even in that case the expected event rates are not very large - and we haven't 
included any efficiencies yet. Optical observations of air showers are typically only possible 
on clear, moonless nights, i.e. some 10% of the time. Once this is included, only the most 
optimistic scenario (a homogeneous distribution of 10 12 GeV particles, and an LSP-nucleon 
scattering cross section near the upper end of the expected range) will give a clear signal in 
this kind of experiment. However, an enlargement of the target area is at least conceivable, 
e.g. by simply launching more satellites, or by using a higher orbit (with correspondingly 
tougher requirements on the optics). 

Table 2: Number of LSP events with E% > 10 9 GeV emerging from at least 5° below the horizon 
expected each year in space-borne fluorescence detectors, for a^ p = a up /10; a ten times smaller 
scattering cross section gives about 8 to 9 times smaller rates. The notation is as in table 1. 



m x [GeV] 


Galactic dist. 


Homogeneous dist. 


2 • 10 12 
2 • 10 16 


2-30 
0.2-0.6 


30-400 
3-10 



4 Summary and Conclusions 

In this contribution I discussed particle physics explanations of the post-GZK cosmic ray 
events. My personal favorite remains the explanation in terms of superheavy late decaying 
particles, which can solve both the energy and propagation problems. The biggest current 
challenge of this class of models seems to be how to explain the paucity of photons as primaries 
in the observed events. A first minnowing of models should be achieved by measuring the flux 
of neutrinos with E u > 10 5 GeV. A decisive test may require measuring (or excluding) the 
predicted flux of LSPs with E^ > 10 9 GeV, which will be difficult, but may be possible using 
space-based experiments. Proving the existence of such mysterious super-heavy but long-lived 
particles would revolutionize our understanding of physics. If these models are still viable after 
the next generation of cosmic ray observatories and neutrino telescopes have presented their 
data, no effort should be spared to search for these extremely energetic cosmic LSPs. 
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